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1. Introduction 

Fractional calculus is a discipline that studies integrals and derivatives of non- 
integer (real or complex) order ^3jil8, 27, 34^. The subject is nowadays very ac- 
tive due to its many applications in mechanics, chemistry, biology, economics, 
and control theory [36] . In 1993, Samko and Ross proposed an interesting gen- 
eralization to fractional calculus |35j (see also Samko's paper |^ of 1995). They 
introduced the study of fractional integration and differentiation when the order 
is not a constant but a function. More precisely, they considered an extension of 
Riemann-Liouville and Fourier definitions [32, .33.,35 . Afterwards, several works 
were dedicated to variable order fractional operators, their applications and in- 
terpretations (see, e.g., [HITIUS])- In particular, Samko's variable order fractional 
calculus turns out to be very useful in mechanics and in the theory of viscous 
flows [7l[9l[T9l[26l[28l[29] . Indeed, many physical processes exhibit fractional-order 
behavior that may vary with time or space |19| . The paper is devoted to the 
study of a variable-order fractional differential equation that characterizes some 
problems in the theory of viscoelasticity. In [5] the authors analyze the dynamics 
and control of a nonlinear variable viscoelasticity oscillator, and two controllers 
are proposed for the variable order differential equations that track an arbitrary 
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reference function. The work [26] investigates the drag force acting on a particle 
due to the oscillatory flow of a viscous fluid. The drag force is determined using 
the variable order fractional calculus, where the order of derivative vary according 
to the dynamics of the flow. In |29) a variable order differential equation for a 
particle in a quiescent viscous liquid is developed. For more on the application of 
variable order fractional operators to the modeling of dynamic systems, we refer 
the reader to the recent review article [28] . 

In this note we develop the fractional calculus of variations via a variable 
order approach. The fractional variational calculus was born in 1996-1997 with 
the works of Riewe in mechanics [301131] , and is now under strong current research 
(see [?ll51[51I^IFllTUllTT]rHHT71[Mf25] and references therein). However, to the best of 
the authors knowledge, results for the variable order case are a rarity and reduce to 
those in [4 . Motivated by the advancements of [31[53] , and in contrast with , we 
consider here fractional problems of the calculus of variations where the Lagrangian 
depends on classical integer order derivatives and both on variable order fractional 
derivatives and integrals. 

The paper is organized as follows. In i|2] a brief review to the variable order 
fractional calculus is given. Our results are then formulated and proved in i]3] we 
show the boundedness of the variable order Riemann-Liouville fractional integral 
in the space Li[a, b] (Theorem l3.ip : integration by parts formulas for variable order 
fractional operators (Theorems 13. 21 and 13. 3p : and a necessary optimality condition 
for our general fundamental problem of the variable order fractional variational 
calculus (Theorem 13. 5p . Two illustrative examples are discussed in ij4j 

2. Preliminaries 

The reference book for fractional analysis and its applications is [34] . Here we recall 
the necessary definitions for the variable order fractional calculus (see, e.g., [19]). 

Definition 2.1 (Left and right Riemann Liouville integrals of variable order). Let 

< a{t, r) < 1 for aU t,T e [a, b] and / G £i [a, b]. Then, 

t 

<^^^'''^f(t)- I n I ,A t~Tr'^'^^^~\f{T)dT it>a) 

J r(a(i,T)) 

a 

is called the left Riemann-Liouville integral of variable fractional order «(•, •), 
while 

b 

J T{a{T,t)) 
t 

denotes the right Riemann-Liouville integral of variable fractional order q;(-, •). 

Definition 2.2 (Left and right Riemann Liouville derivatives of variable order). 

Let < a(i,T) < 1 for all i,T e [a, 6]. If a4""^'^V e AC[a,b], then the left 
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Riemann-Liouville derivative of variable fractional order «(•,-) is defined by 

a 

while the right Riemann-Liouville derivative of variable fractional order a(-, ■) is 
defined for functions / such that til AC[a, b] by 

t 

Definition 2.3 (Left and right Caputo derivatives of variable fractional order). Let 

< a{t,T) < 1 for all t,T G [a,b]. If / G v4C[a, 6], then the left Caputo derivative 
of variable fractional order a(-, •) is defined by 

t 

a 

while the right Caputo derivative of variable fractional order a(-, •) is given by 

b 

t 

The following result will be useful in the proof of Theorems 13.11 and 13.21 
Theorem 2.4 (cf. [12]). Letx G [0, 1]. The Gamma function satisfies the inequalities 

x + l - ^ - a; + 2 ^ ' 



3. Main results 

In ij3.1l we prove that the variable fractional order Riemann-Liouville integral 
aJ^^ ' ^ is bounded on the space Li[a, b] (Theorem l3.ip . In §3. 21 we obtain a formula 
of integration by parts for Riemann-Liouville integrals of variable order (Theo- 
rem [321) a formula of integration by parts for derivatives of variable fractional 
order (Theorem 13. 3p . Finally, in ij3.3l we use the obtained formulas of integration 
by parts to derive a necessary optimality condition for a general problem of the 
calculus of variations involving variable order fractional operators (Theorem [33]). 

3.1. Boundedness 

The following result allow us to consider problems of the calculus of variations with 
a Lagrangian depending on left Riemann-Liouville integrals of variable order. 
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Theorem 3.1. Let i < a{t, t) < 1 for all t,T G [a, b] and a certain n e N greater 

or equal than two. The Riemann-Liouville integral ait'^ ' : Li[a,b] Li[a,b] of 
variable fractional order •) is a linear and bounded operator. 



Proof. The operator is obviously linear. Let i < a{t,T) < 1, / 6 ii[a, &], and 
define 

I if r > i 



F{T,t) 



for all (r, i) e A = [a, fe] x [a,b]. Since i < a{t,T) < 1, then 

1. for T + 1 < t we have ln(t - r) > and (t - t)"(*^^)-i < 1; 

2. for T < t < T + 1 wc have ln(< - t) < and {t - t)°'(*'^)-'^ < {t - r) tt-i. 

Therefore, 



a \ a 



F{T,t)dt\dT^ / |/(r) 



r(«(t,r)) 



a(t,T)~l 



= / 



_ _^a(t,T)-l 



< / I/Ml 



dt 



T + l 



(t-r) 



r(a(i,r)) 



dt dr 



dt dr 



r(a(t,T)) 

Moreover, by inequality (|2.ip . one has 





1 




dt+ [ 




r(a(t,T)) 


Jt+1 



l/WI 



< / I/Ml 



r+l 



r(a(i,r)) 





1 




dt+ [ 

Jt+1 


j 


r(«(t,r)) 



a^{t,T) + l 



{t-T)^^^dt + 



T+l 



a ^{t,T) +a{t,T) 
a2(t,T) + 1 



dt dr 



< l/(r)| 



/ {t -T)i-^dt + b~T -l\ ch = \f{T)\{n + b-T -l)d7 

\.J T / J a 

< + 11/11 <oo. 

It follows from Fubini's theorem that F is integrable in the square A and 



1 



< 



a r(a(i,T)) 
6 / 1 



\a(t,T)-l 



/(r)dr 



a \>/ a 



r{a{t,T)) 



dt 

dr] dt 



(^J^ F{T,t)dT^ dt 



<{n + b-a)\\f\\. 
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Therefore, ^ ' "* : Li[a, b] — > Li[a, b] and a^t 



»(■,■) 



< n + b — a. 



□ 



3.2. Integration by parts formulas 

The integration by parts formulas we now obtain have an important role in the 
proof of the generalized Euler-Lagrange equation (I3.5p . We note that in Theo- 
rem [321 the left-hand side of p. II) involves a left integral of variable order while 
on the right-hand side it appears a right integral. 

Theorem 3.2. Let ^ < a{t, t) < 1 for all t,T e [a, b] and a certain n G N greater 
or equal than two, and f,g E C ([a, &];M). Then, 



\{t)ai^^-y{t)dt= f f[t)tr,^-\{t)dt. 



(3.1) 



Proof. Define 

F{T,t) 



i{t.T)-l 



9{t)f{r) 







if T < i; 
if T > i 



for all (r, t) G A = [a,b] x [a,b]. Since / and g are continuous functions on [a, 6], 
they are bounded on [a,b], i.e., there exist Ci,C2 > such that \g{t)\ < Ci and 
\f{t)\ <C2,te [a,b]. Therefore, 



f I f F{T,t)dT\ dt^ f 



T{a{t,T)) 



it 



_ _^a(t,r)-l 



dr dt 



a \J a 
b 



-it-r) 



Q(t,T)-l 



dr dt 



T{a{t,T)) 

^'^^r(/'rwb))<'-^'°"'"""*)''- 

Since ^ < a{t,T) < 1, then 

1. for 1 < t - T we have ln(< - r) > and {t - r)"(*^^)-i < 1; 

2. for 1 > t - T we have ln(< - t) < and {t - r)"(*'^)-i < (t - r) 

Therefore, 



C1C2 



1 



'a r(a(i,T)) 

.6 / rt~l 

<CiC2 



(t_r)"(*'^)-idr ) di 
1 



T{a{t,T)) 



dr 



1 



t-i 



r(a(t,r)) 



(< - T)""^dr ) dt. 
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Moreover, by (12.11) . one has 



C1C2 



t-i 



1 



' a \'J a 
< C1C2 



r(a(i,r)) 

'-1 ^,2 



dT 



1 



t-i 



r(a(t,r)) 



(t-r)-"VT dt 



a^{t, r) + a{t, r) 



dT 



a^{t, r) + a{t, r) 
a^{t,T) + 1 



= CiC2ib-a) 



I { I dT+ f {t-T)i-^dT] dt 
Ja \Ja Jt-1 J 



dT 

b + a 



1 + n — a ] < 00 



Hence, one can use the Fubini theorem to change the order of integration: 



9{t)f{r) 



1 



T{a{t,T)) 



it - r)"(*--)-ldi ) dT = f f{T)rt-^-\{T)dT. 



□ 

In our second formula p.2p of fractional integration by parts, the left-hand 
side contains a left Caputo derivative of variable fractional order q;(-, •), while on 
the right-hand side it appears a right Riemann-Liouville integral of variable order 
1 — a(-, •) and a right Riemann-Liouville derivative of variable order q;(-, •). 



Theorem 3.3. Let Q < a{t,T) < 1 — — for all t,T E [a, b] and a certain n G N greater 



or equal than two. If f € ([a,6];M), g G C{[a,b];R), and jI "^'''''g € AC[a,b], 



g{tf,Df-'-^f{t)dt^ f{t)al'^^-'-^g{t) + f{t)tDf-'-^git)dt. (3.2) 



then 



Proof By Definition O it follows that '^Df''''> f(t) = al\ ''^ ft f ■ Applying 
Theorem 13.21 and integration by parts for classical (integer) derivatives, we obtain 



g{t)^Dt^-^-^fit)dt: 



9{t) 



m,fr^-'-^git) 



b 



g{t)dt 



d 



f{t)i-A-''^-'-^9{t)dt 



= f{t)tfr'^-'\{t) + / f{t),Di'^-'-\{t)dt. 



□ 
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3.3. A fundamental vjiriational problem of vstriable fractionsJ order 

We consider the problem of extremizing (minimizing or maximizing) a functional 

b 

J[yi-)] = / ^ y{t),y'it),^Df-'-\it), (-^2/(0) dt (3.3) 

a 

subject to boundary conditions 

y{a) = Va, y{b) = Vb, (3.4) 

where a,/3 : [a,b] x [a,b] — > M are given functions taking values in (0,1 — i) 
and 1), respectively, with n G N greater or equal than two. For simplicity of 
notation, we introduce the operator {•,-,•} defined by 

{j/,a,/3} {t) = (i,y(t),y'(t),^Dr(-'\(0,aif<-^\(i)) • 

We assume that F e M*;R); t i-^ d^F {y,a, (3} (t) is continuous and 

has absolutely continuous integral til ' ' and continuous derivative t-D^^ ' ^; 1 1— ^ 
dsF {y, a, /3} (t) is continuous and has continuous variable order fractional integral 
tjf*' ' and t H> dsF {y, a, /?} {t) has continuous usual derivative ^. 

Definition 3.4. A continuously differentiable function y ^ ([a, 6];R) is said to 

be admissible for the variational problem p.3P " P.4p . if ^D"^ ' ^y and ajf ^ ' ''y ex- 
ist and are continuous on the interval [a,b], and y satisfies the given boundary 
conditions (|3.4p . 

Theorem 3.5. Let y be a solution to problem (|3.3p - p.4p . Then, y satisfies the 
generalized Euler-Lagrange equation 

d^F {y, a, /?} (t) - ^^3^ {y, a, /?} {t) + -^SsF {y, a, /?} (t) 

+ tD^('-)94F{y,a,/3}(i)=0 (3.5) 

/or all t G [a, 6]. 

Proof. Suppose that y is an extremizer of . Consider the value of ^7 at a nearby 
function y{t) — y{t)+eri{t), where e S M is a small parameter and rj G {[a, b];M.) 
is an arbitrary function satisfying 77(a) = 77(6) — and such that a^^t^ ' ^vi^) s-^d 
' ^y{t) are continuous. Let 

J{e) = J[y{-)] = I F{y,a,l3}{t)dt. 
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A necessary condition for y to be an extremizer is given by 



dJ 

de 



6=0 



^ lid^F {y, a, /?} {t) ■ rjit) + d^F {y, a, /?} {t)j^v{t) 



+ d^F {y, a, /?} {tf,Df-'\{t) + d,F {y, a, /?} (t) • at''' ''^ v{t)j dt = 0. (3.6) 

Using the classical and the generalized fractional integration by parts formulas of 
Theorems 13.21 and 13. 3[ we obtain 

h b 

J diF'^Df-'-^r,dt=r^tll-''^-'-^diF\\ J VtOf''^ O.Fdt, 



and 

6 6 

J d^Fj,^'\dt = J rjtC'-^d^Fdt. 

a a 

Because 77(a) — rj{h) — 0, ()3.6p simplifies to 

6 / . 
77(0 I d^F {y, a, /?} {t) - -d^F {y, a, /3} (i) + tDl^'^ '^dAF {y, a, /?} {t) 

+ ti^''^d^F{y,a,P} {t)^dt = Q. 

One obtains (|3.5p by the fundamental lemma of the calculus of variations (see, 
e.g., [37]). □ 

4. Illustrative examples 

Let /3(t, r) = j3{t) be a function depending only on variable t, ^ < (3{t) < 1 for all 
t £ [a, b] and a certain 77 G N greater or equal than two, and 7 > — 1. We make 
use of the identity 

' ^ ^ - r(7 + /3(t) + i) ^^-^^ 

that one can find in the Samko and Ross paper |35) . 
Example 1. Let J7 be the functional defined by 



= / V^' Sr(3,r.')''«') (-'^'''''»)'-°'^'''''' 



dt 
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with endpoint conditions y{a) = and y{b) = (6 — a)^. If y is an extremizer for 
then the necessary optimality condition of Theorem 13.51 gives 



7/3(-) 



( r(/3(t)+3) f^■),,u\_^ ^ 

r(3)(t-a)2+/3(*) a^t yV-l ^ 



By identity (|4.ip . function 



= 0. (4.2) 



y(t) = (i - a)^ (4.3) 
is a solution to the variable order fractional differential equation 

^(.) _ r(3)(t-«)^+/^(*) 
^ r(/?(t) + 3) ■ 

Therefore, ()4.3p is a solution to the Euler-Lagrange equation (|4.2p and an extremal. 

In the next example a(i,T) is a function taking values in the set (0, 1 — -i). 
Example 2. Consider the following problem: 

b 2 

J[y{-)]- J {^Df-^y{t)y + (aif^\(t) - r('/3ft)"+ 1) ) '^^^^i^' 

a 

?/(«) = e, y(6) = e, 

for a given real ^. Because J'[2/(-)] > for any function y and J'[y{-)] — for the 
admissible function y — ^ (use relation (|4.1I) for 7 = 0, the linearity of operator 
alt^ \ and the definition of left Caputo derivative of a variable fractional order), 
we conclude that y is the global minimizer to the problem. It is straightforward to 
check that y satisfies our variable order fractional Euler-Lagrange equation p.Sp . 
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